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Abstract. We present the results of a high spectral resolution (A/AA = 49000) study of the circumstellar (CS) gas around the 
intermediate mass, pre-main sequence star UX Ori. The results are based on a set of 10 echelle spectra covering the spectral 
range 3800 - 5900 A, monitoring the star on time scales of months, days and hours. A large number of transient blueshifted 
and redshifted absorption features are detected in the Balmer and in many metallic lines. A multigaussian fit is applied to 
determine for each transient absorption the velocity, v, dispersion velocity, Av, and the parameter R, which provides a measure 
of the absorption strength of the CS gas. The time evolution of those parameters is presented and discussed. A comparison of 
intensity ratios among the transient absorptions suggests a solar-like composition of the CS gas. This confirms previous results 
and excludes a very metal-rich environment as the cause of the transient features in UX Ori. The features can be grouped by 
their similar velocities into 24 groups, of which 17 are redshifted and 7 blueshifted. An analysis of the velocity of the groups 
allows us to identify them as signatures of the dynamical evolution of 7 clumps of gas, of which 4 represent accretion events 
and 3 outflow events. Most of the events decelerate at a rate of tenths of m s~ 2 , while 2 events accelerate at approximately the 
same rate; one event is seen experiencing both an acceleration and a deceleration phase and lasts for a period of few days. 
This time scale seems to be the typical duration of outflowing and infalling events in UX Ori. The dispersion velocity and the 
relative aborption strength of the features do not show drastic changes during the lifetime of the events, which suggests they 
are gaseous blobs preserving their geometrical and physical identity. These data are a very useful tool for constraining and 
validating theoretical models of the chemical and physical conditions of CS gas around young stars; in particular, we suggest 
that the simultaneous presence of infalling and outflowing gas should be investigated in the context of detailed magnetospheric 
accretion models, similar to those proposed for the lower mass T Tauri stars. 

Key words. Stars: formation - Stars: pre-main sequence - Circumstellar matter - Accretion, accretion disks - Lines: profiles 
- Stars: individual: UX Ori 



1. Introduction 
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is nowadays accepted that planets form from CS disks via 
the formation of planetesimals (Beckwith et al. 2000). Several 
lines of evidence suggest that the young main sequence A5V 



Table 1. EXPORT UES/WHT UX Ori observing log. The 
Julian date (-2450000) of each exposure is given in column 1 . 
Columns 2 to 5 give simultaneous V, H, K and %P V photopo- 



star p Pic (20 Myr, Barrado y Navascues et al. 1999) hosts larimetric data, when available, taken from Oudmaijer et al 



planetesimals inside its large CS disk. The main argument is 
the presence of transient Redshifted Absorption Components 
(RACs) in high resolution spectra of strong metallic lines, like 
Can K 3934 A. These spectroscopic events have been inter- 
preted as being caused by the evaporation of comet-like highly 
hydrogen-depleted bodies. The interpretation is known as the 



Falling Evaporating Bodies (FEB) scenario (Lagrange et al 



200C, and references therein). However, the presence or ab- 
sence of planetesimals during the pre-main sequence (PMS) 
phase is a controversial observational to pic. The time scale for 
the formation of planetesimals (~10 4 yr, Beckwith et al. 200C ) 
is shorter than the duration of the PMS phase (~ 1—10 Myr), 
which suggests that they should exist during PMS stellar evo- 
lution. 

UX Ori-like PMS objects (UXORs) are characterized by a 
peculiar photo-polarimetric variability, which has been inter- 
preted as the signa ture of massive, almost edge-on, CS disks 
( Grinin et al. 1991 ). Most UXORs have A spectral types and 
therefore are the PMS evolutionary precursors of f3 Pic. UXORs 
have been reported to show RACs (Grinin et al. 1994; de 
Winter et al. 1999); they also show Blueshifted Absorption 
Components (BACs) in their spectra. In this paper, the acronym 
TAC (Transient Absorption Component) will be used to denote 
both RACs and BACs. In analogy to p Pic, RACs observed in 
UXORs have been interpreted in terms of the FEB scenario 
(Grady et al. 2000). However, this interpretation has recently 



been questioned by Natta et al.| (2000), who used the spectra 
obtained by prinin et~aT (2001) to analyze the dynamics and 
chemical composition of a very strong, redshifted event in UX 
Ori itself, an A4 IV star (|Mora et al.||2001|) ~ 2 x 10 6 year 



old ( Natta et al. 1999| ). Gas accretion from a CS disk was sug 



gested in Matta et al. [( |2000| ) as an alternative to the FEB sce- 
nario to explain the observed RACs in UX Ori. In addition, 



Beust et al 



( 2001 ) have found that the FEB hypothesis can- 
not produce detectable transient absorptions in typical HAe CS 
conditions, unless the stars are relatively old (> 10 7 yr). 

A detailed observational study of the kinematics and chem- 
istry of TACs, e. g. by means of echelle spectra which simulta- 
neously record many metallic and hydrogen lines, can discrim- 
inate between the two scenarios. For instance, in a FEB event 
large metallicities are expected, while gas accreted from a CS 
disk would have approximately solar abundances. A strong ob- 
servational requirement is set by the time scale of monitoring of 
the TACs. In [Natta et aL ( 2000 ) spectra were taken 3 days apart. 
Since UX Ori is a highly variable star, there is some ambiguity 
in identifying transient spectral components of different veloc- 
ities detected over this time interval as having the same physi- 
cal origin. A better time resolution is needed in order to ensure 
that the TACs observed at different velocities are due to the 
dynamical evolution of the same gas. The EXPORT collabora- 



tion (Eiroa et al. 2000) obtained high resolution echelle spec- 



tra of a large sample of PMS stars (Mora et al. 2001). About 
10 PMS stars showed TACs which were intensively monitored 
(At < 1 day). The study of the kinematics and chemistry of 
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these events provides important tools for identifying their ori- 
gin, or at least to put severe constraints on it. 

This paper presents an analysis of the TACs observed in 
UX Ori by EXPORT and shows that the results are not com- 
patible with a FEB scenario. The layout of the paper is as fol- 
lows: Section 2 presents a brief description of the EXPORT 
observations. Section 3 presents the procedures followed in the 
analysis of the spectra. Section 4 presents the kinematic and 
chemical results of the detected TACs. Section 5 gives a dis- 
cussion of the dynamics and nature of the gas. In section 6, we 
present our conclusions. 

2. Observations 

High resolution spectra of UX Ori were taken in October 
1998 and January 1999 using the Utrecht Echelle Spectrograph 
(UES) at the 4.2m WHT (La Palma Observatory). We collected 
10 echelle spectra in the wavelength range 3800 - 5900 A, 
with resolution /1/A/l = 49000 (6 km/s); exposure times range 
between 20 and 30 minutes. The observing log is given in 
Table [j], in which the long-term (months) and short-term (days, 
hours) monitoring is evident. Standard MIDAS and IRAF pro- 
cedu res ha ve been used for the spectroscopic reduction (Mora 
et al. 2001 ). Final typical signal to noise ratio (SNR) values are 
-150. 

EXPORT also carried out simultaneous intermediate res- 
olution spectroscopy, optical photo-polarimetry and near-IR 
photometry of UX Ori during the nights when the UES spec- 
tra were taken ( [Vlora et aT 2001; Oudmaijer et al. 2001; 



Eiroa 



et al. |2001| ). Table |l] also gives photo-polarimetric data taken 



simultaneously with the UES spectra; they show that UX Ori 
was always in its bright state with little variation, but significant 
(|Oudmaijer et al.||200l[). 



3. Analysis of the spectra 

A detailed analysis of TACs requires high SNR (> 25) and high 
spectral resolution (< 10 km/s) in order to resolve the kine- 
matic components observed simultaneously in different lines, 
as is shown in Fig. 1. This figure presents the Nai D lines in 
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Fig. 1. Comparison between intermediate (top, J.D. 
2451210.3098, 10 min. exposure time) and high (bottom, 
J.D. 2451210.3317, 30 min. exposure time) resolution spectra. 
The Na i D spectral region is shown. A sharp IS zero velocity 
component and a much broader RAC (indicated in the figure) 
can be seen in both lines of the doublet in the high resolution 
spectrum. They are also detected in the intermediate resolution 
spectrum, but the big difference in the width of the compo- 
nents is not clear. For the sake of clarity an artificial vertical 
displacement of 0.2 units has been applied to the intermediate 
resolution spectrum. 



one of the UES spectra together with the simultaneous inter- 
mediate resolution spectrum (R ~ 6000) in the same lines. The 
UES spectrum clearly shows a sharp, IS, narrow zero velocity 
component and a RAC in both lines of the doublet while both 
components, though visible, cannot be cleanly separated in the 
lower resolution spectrum. This is why we restrict our further 
analysis of the absorptions to the UES spectra. 

All the UX Ori UES spectra show CS spectral features 
in a variety of lines. The features are seen in absorption and 
can be either blueshifted or redshifted; some underlying emis- 
sion is also present (see below). Their profiles are complex and 
blended components are directly seen in many cases. Thus, the 
analysis of the CS contribution to the observed spectra requires 
a careful subtraction of the stellar photospheric spectrum and 
a method to characterize the blended kinematical components. 
In this section we descr ibe the procedure we have followed, 
which is the one used in Natta et al. [ ( |2000b , although we have 
improved it by considering all detected TACs simultaneously. 



3.1. Subtraction of the UX Ori photospheric spectrum 

Firstly we estimate the radial velocity of the star. The helio- 
centric correction of each observed spectrum is computed us- 
ing MIDAS. Then, the radial velocity of the star is estimated 
by fitting the Na I D sharp IS absorption. This absorption is 
stable in velocity (velocity differences in our spectra are less 
than 1.0 km/s) and is related to the radial velocity of the star 
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Fig. 2. Comparison of the observed average spectrum of UX 
Ori (solid line) to the synthetic spectrum with stellar parame- 
ters T e ff = 9250 K, log g = 4.0 and v sin i — 215 km/s (dashed 
line). The spectral region shown includes the blends at 4172 - 
4179 A (mainly Fe n and Ti n) and 4203 A (Fe i). 



is obtained, in good agreement with the 18 km/s estimate by 
|Grinin et al.| ( |l994| )- 

Kurucz (1993) model atmospheres assuming solar metal- 
licity and turbulence velocity of 2 km/s have been used to 
synthesize the photospheric spectra, which are later used for 
a comparison with the observed ones. Atomic line data have 



been obtained from the VALD online database (Kupkaetal 



1999 ), and the codes SYNTHE ([Kurucz |1979| ; pfeffery et al 
1996] ) and SYNSPEC ( |Hubeny et al.||l995| ) are used to synthe 



size the metallic lines and Balmer hydrogen lines, respectively. 
In this way we compute a large grid of synthetic spectra with 
the effective temperature, gravity and rotation velocity as free 
variable parameters. 

The best set of free parameters is estimated by comparing a 
number of lines among the observed spectra and the synthetic 
ones. This is not an easy task since most UX Ori lines are vari- 
able to some extent, and the choice of pure photospheric lines 
is not trivial. Our choice is to consider a number of weak lines 
that show a minimum degree of variability (less that 1%) in 
the UES spectra. The blends located at 4172 - 4179 A (mainly 
Fe ii and Ti n) and 4203 A (Fe i) have been used to estimate 
T e ff (effective temperature) and log g (logarithm of the surface 
gravity), since they are sensitive to changes in T e ff and log g 
for early type A stars (Gray & Garrison 1987; Gray & Garrison 



1989). The synthetic spectra have also been broadened allow- 



( Finkenzeller & Jankovics 1984). A value of 18.3 + 1.0 km/s 



ing for stellar rotation. The best set of free parameters we find 
is T e ff = 9250 K, log g = 4.0, v sin i — 215 km/s. Fig. || shows 
the comparison between the synthetic and the average observed 
spectrum in the above-mentioned lines; the agreement is excel- 
lent. The stellar parameters of the synthetic spectrum are in 
very good agreement wit h the spectral typ e and rotational ve- 
locity derived elsewhere ( Mora et al. 2001 ). 

The next step consists of the subtraction of the best syn- 
thetic photospheric spectrum from each observed spectrum. 
This has been carried out by fitting the continuum with a lin- 
ear law between two points selected at both sides of those lines 
showing TACs. As an example Fig. shows the Ca n K line re- 
sults, where the synthetic spectrum and the J.D. 2451209.4204 
UX Ori observed one are plotted together with the residual af- 
ter the subtraction. The residual corresponds to the CS contri- 
bution to the observed spectra. As a more convenient way to 
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Fig. 3. Subtraction of the photospheric spectrum. The spec- 
tral range corresponds to the Can K line. Top: observed 
(solid line) and synthetic (dotted line) spectra are displayed. 
Middle: Residual due to the CS contribution (observed - pho- 
tospheric). Bottom: Normalized residual absorption, R = (1 - 
observed / photospheric). Vertical tick marks indicate the 4 ab- 
sorption components of the Ca n K which can be identified in 
this spectrum, including the narrow IS absorption. In the bot- 
tom we show the 4 gaussians (dashed lines), whose sum (dotted 
line) excellently reproduces the R profile (solid line). 



represent this contribution we define the R parameter (normal- 
ized residual absorption): 

R - I - F bs/F syn 

R — means no CS absorption, and R = 1 denotes com- 
plete stellar light occultation, i.e. R quantifies the CS absorp- 
tion strength and could be directly compared to the models in 
Natta et al.| ( ]2000| ), R is also plotted in Fig. |j 



3.2. R multicomponent gaussian fitting 

The CS contribution presents a very complex line profile 
caused by the blending of different transient features. As an ex- 
ample, Fig. |] shows the components of the Can K line which 
are clearly present in that particular observation. For the sake 
of analysing the profile we assume that it is caused by gas 



with different kinematics and that each component can be rep- 
resented by means of a gaussian function. This choice pro- 
vides a direct interpretation of the fit parameters in terms of 
a physical characterization of the gaseous transient absorption 
components: the gaussian center gives the central velocity (V) 
of each component, the gaussian FWHM (Full Width at Half 
Maximum) provides the velocity dispersion (Av) and finally 
the gaussian peak gives the value of the residual absorption 
R. Since the feaures are blended a multigaussian fit is needed, 
which has been carried out by means of the IRAF ngaussf it 
routine. Fig. ^ shows the results of the multigaussian fit for 
the Can K line. 4 gaussians, i.e. 4 kinematical components, 
each with their corresponding parameters, provide an excellent 
agreement with the residual CS contribution profile. 

This analysis is applied to the lines with the strongest vari- 
ability, which are good tracers of the CS gas properties: H/3 
4861 A, Hy 4340 A, U6 4102 A, He 3970 A, H£ 3889 A, Ca n K 
3934 A, Ca n H 3968 A, Na i D2 5890 A and Na i D 1 5896 A. 
All Balmer lines with good SNR present transient absorption 
features (even the noisy Hk 3750 A, 10 th line in the Balmer se- 
ries). They are also observed in man y other lines (e.g. H e i, Fe i 
and Fen), as already pointed out by Grinin et al. (2001), but a 
corresponding discussion is deferred to a future work. 

Underlying line emission is often significant in the Balmer 
lines, particularly in H/3, Hy and H6; in fact, Ha is seen in emis- 
sion (equivalent width ~ 10 A) in the simultaneous INT spec- 
tra. Ca ii K and H also show little emission. In these cases R is 
underestimated and negative values can artificially arise when 
the synthetic photospheric spectrum is subtracted. This is es- 
pecially bad for faint absorption components near the emission 
peaks since the latter appear as minima in the R plot. In order 
to avoid contamination by these emission components, we de- 
fine a new "zero level", estimated by a linear fit to the minima 
observed in the R plot. Fig. |] illustrates this approach. The con- 
sistency of the results obtained with different lines confirms its 
validity. Further improvements should be based on theoretical 
modelling or higher order fitting of the emission. 

4. Results 

Table || gives the radial velocity shift v, the velocity dispersion 
Av and the absorption strength R of each identified transient ab- 
sorption, computed according to the procedure described in the 
previous section. Uncertainties affecting the values of Table || 
are difficult to quantify, especially in the case of blended com- 
ponents. We are confident, however, that the values are a good 
representation of the gas properties in the case of sharp, iso- 
lated features. The results allow us to identify trends, based on 
the similar and consistent behaviour of many different lines. 

Absorption components at different radial velocities - ei- 
ther redshifted or blueshifted - are found in the Balmer and 
metallic lines. Can and Nai components with the stellar radial 
velocity are also detected. Features in different lines with sim- 
ilar radial velocities are detected within each spectrum; thus, 
it is reasonable to assume that they form in approximately the 
same region. The absorptions with similar velocities appearing 
in different lines are called a TAC, which can be characterized 
by the average of the radial velocity of the individual lines. A 
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Table 2. Identified transient absorption components in the Balmer and Can and Nai lines of UX Ori. JD is the Julian Date 
(-2450000). Numbers in column 3 represent the event assigned to the particular absorptions (see text section 4). A "0" in that 
column corresponds to the narrow IS absorptions, while a "-" means that the absorption is not associated with a particular event. 
Columns 4 to 6 give: v, the radial velocity of the transient absorption, Av, the FWHM and R, the normalized residual absorption. 
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4 


138 


287 


0.66 


Nai D2 


1113.6034 


2 


61 


38 


0.09 


Hy 


1210.3317 


4 


138 


255 


0.54 


H/J 


1113.6034 


3 


-51 


42 


0.07 


Ha 


1210.3317 


4 


142 


270 


0.43 


Hy 


1113.6034 


3 


-30 


69 


0.07 


He 


1210.3317 


4 


110 


277 


0.64 


CauK 


1113.6034 


3 


-56 


45 


0.06 


Hf 


1210.3317 


4 


151 


265 


0.31 


CanK 


1113.6034 





2 


17 


0.30 


CanK 


1210.3317 


4 


124 


265 


0.54 


CanH 


1113.6034 





5 


15 


0.26 


Nai D2 


1210.3317 


4 


115 


137 


0.06 


Nai D2 


1113.6034 





5 


12 


0.68 


Nai Dl 


1210.3317 


4 


133 


150 


0.06 


Nai Dl 


1113.6034 





5 


12 


0.53 


H/J 


1210.3317 


7 


3 


46 


0.19 


H/J 


1113.7194 


1 


114 


188 


0.36 


Hy 


1210.3317 


7 


9 


39 


0.10 


Hy 


1113.7194 


1 


90 


169 


0.36 


CauK 


1210.3317 





-3 


19 


0.20 


m 


1113.7194 


1 


182 


89 


0.11 


CanH 


1210.3317 





1 


8 


0.16 


He 


1113.7194 


1 


143 


121 


0.20 


Nai D2 


1210.3317 





1 


12 


0.57 


Hi 


1113.7194 


1 


165 


84 


0.09 


Nai Dl 


1210.3317 





2 


13 


0.46 


CauK 


1113.7194 


1 


84 


205 


0.35 


H/J 


1210.3568 


4 


132 


287 


0.64 


H/J 


1113.7194 


2 


52 


75 


0.62 


Hy 


1210.3568 


4 


132 


260 


0.53 


Hy 


1113.7194 


2 


57 


58 


0.44 


Ha 


1210.3568 


4 


128 


238 


0.43 


H<5 


1113.7194 


2 


59 


121 


0.56 


He 


1210.3568 


4 


107 


280 


0.63 


He 


1113.7194 


2 


50 


96 


0.51 


Hf 


1210.3568 


4 


140 


249 


0.29 


h<t 


1113.7194 


2 


57 


108 


0.38 


CanK 


1210.3568 


4 


119 


296 


0.52 


CauK 


1113.7194 


2 


56 


50 


0.35 


Nai D2 


1210.3568 


4 


113 


109 


0.04 


CanH 


1113.7194 


2 


60 


60 


0.60 


Nai Dl 


1210.3568 


4 


142 


127 


0.04 


Nai D2 


1113.7194 


2 


61 


44 


0.10 


H/J 


1210.3568 


7 


6 


41 


0.18 


Nai Dl 


1113.7194 


2 


60 


28 


0.06 


Hy 


1210.3568 


7 


17 


37 


0.08 


H/J 


1113.7194 


3 


-44 


72 


0.13 


CauK 


1210.3568 








15 


0.14 


Hy 


1113.7194 


3 


-22 


99 


0.13 


CanH 


1210.3568 





-2 


12 


0.14 


CanK 


1113.7194 


3 


-58 


52 


0.08 


Nai D2 


1210.3568 





2 


13 


0.60 


CanK 


1113.7194 





6 


19 


0.31 


Nai Dl 


1210.3568 





3 


13 


0.48 


CanH 


1113.7194 





5 


19 


0.20 


H/J 


1210.4237 


4 


134 


240 


0.66 


Nai D2 


1113.7194 





5 


13 


0.64 


Hy 


1210.4237 


4 


122 


242 


0.55 


Nai Dl 


1113.7194 





5 


12 


0.49 


Ha 


1210.4237 


4 


125 


229 


0.44 


H/J 


1207.5268 


4 


78 


221 


0.26 


He 


1210.4237 


4 


89 


253 


0.64 


Hy 


1207.5268 


4 


48 


211 


0.29 


Hf 


1210.4237 


4 


126 


230 


0.30 


US 


1207.5268 


4 


52 


223 


0.24 


CauK 


1210.4237 


4 


117 


272 


0.51 


He 


1207.5268 


4 


68 


207 


0.17 


Nai D2 


1210.4237 


4 


120 


99 


0.05 


Hf 


1207.5268 


4 


83 


128 


0.16 


H/J 


1210.4237 


7 


10 


53 


0.22 


CanK 


1207.5268 


4 


51 


269 


0.32 


Hy 


1210.4237 


7 


20 


26 


0.07 


H/J 


1207.5268 


5 


-60 


88 


0.61 


CauK 


1210.4237 





-1 


18 


0.15 


Hy 


1207.5268 


5 


-53 


90 


0.46 


CanH 


1210.4237 





-2 


12 


0.13 


US 


1207.5268 


5 


-49 


82 


0.39 


Nai D2 


1210.4237 





2 


13 


0.61 


He 


1207.5268 


5 


-50 


90 


0.19 


Nai Dl 


1210.4237 





3 


13 


0.50 


Hf 


1207.5268 


5 


-43 


92 


0.30 


H/J 


1210.5156 


4 


124 


217 


0.62 


CauK 


1207.5268 


5 


-59 


70 


0.32 


Hy 


1210.5156 


4 


117 


218 


0.51 


CanH 


1207.5268 


5 


-56 


63 


0.38 


Ha 


1210.5156 


4 


111 


223 


0.40 


Nai D2 


1207.5268 


5 


-54 


22 


0.10 


He 


1210.5156 


4 


76 


262 


0.60 


CauK 


1207.5268 





6 


33 


0.17 


Hf 


1210.5156 


4 


111 


237 


0.25 


CanH 


1207.5268 





1 


21 


0.15 


CauK 


1210.5156 


4 


111 


265 


0.44 


Nai D2 


1207.5268 





3 


13 


0.59 


Nai D2 


1210.5156 


4 


132 


111 


0.04 


Nai Dl 


1207.5268 





3 


12 


0.49 


H/J 


1210.5156 


7 


10 


53 


0.23 


H/J 


1207.5268 




20 


62 


0.20 


Hy 


1210.5156 


7 


13 


49 


0.15 


Hy 


1208.5072 


4 


89 


285 


0.20 


CauK 


1210.5156 


7 


22 


17 


0.11 


HS 


1208.5072 


4 


84 


268 


0.18 


CanH 


1210.5156 


7 


17 


15 


0.07 


He 


1208.5072 


4 


109 


294 


0.14 


CanK 


1210.5156 





-1 


22 


0.20 


Hf 


1208.5072 


4 


91 


331 


0.15 


CanH 


1210.5156 





-3 


15 


0.20 


CauK 


1208.5072 


4 


77 


279 


0.25 


Nai D2 


1210.5156 





3 


13 


0.61 


H/J 


1208.5072 


5 


-35 


107 


0.60 


Nai Dl 


1210.5156 





3 


13 


0.50 
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Fig. 4. Underlying emission in the H/3 absorption line. Top: 
Normalized observed spectrum (solid line) and photospheric 
model (dotted line). The line profile is clearly modulated by 
an emission contribution at =s ±300 km/s in the line wings. 
Middle: Subtracted spectrum (observed - synthetic). Bottom: 
R curve. The vertical mark points to an absorption component 
which would be undetectable if the continuum level were not 
redefined. This has been done by means of a linear fit to the 
peaks of the emission, shown as throughs in the R plot. This fit 
(solid line at the bottom) is the new zero level for the gaussians. 
3 gaussians (dashed lines) whose sum (dotted line) provides the 
optimal fit to the R curve are also plotted at the bottom. 



total of 24 TACs can be identified in our spectra - 17 RACs 
and 7 BACs. Fig. || plots the average radial velocity of each 
detected TAC as a function of the observing time. Error bars 
show the rms error and the number of individual lines used to 
estimate the average velocity is indicated. Fractional numbers 
come from the fact that a weight of 1/2 is assigned to those ab- 
sorption line features that are less certain (the blended lines He 
and Ca n H, the weak line H£ and the Na i doublet affected by 
telluric lines), and weight 1 to the rest. Fig. ^ clearly shows that 
more than one TAC is present in each observed spectrum and 
that there is a radial velocity shift when comparing TACs in dif- 
ferent spectra. The shift does not seem to be accidental; there 
appears to be a systematic temporal evolution of the TACs from 



one spectrum to the following one when their radial velocities 
are analysed. This is clearly evident in the behaviour of the 
TACs observed in the four spectra taken during the night J.D. 
2451210.5 with a time interval of around 4 hours from the first 
to the last one. Thus, the data reveal groups of TACs represent- 
ing a dynamical evolution of the gas with which they are associ- 
ated. We call each of these groups an "event". In total, 7 events 
are identified, of which 4 are redshifted and 3 blueshifted. The 
detailed temporal evolution of v, Av and R for each individual 
line absorption component of 2 of these events (1 blueshifted 
and 1 redshifted) is shown in Fig.^. 

4.1. Kinematics 

In this section we describe the kinematic behaviour of the 
circumstellar gas as revealed by our observations. Since they 
cover only a limited amount of time, it is difficult to estimate if 
the trends we have observed are always present, or if they are 
specific to our observations, and a different choice of times, or a 
longer coverage of the star, would give a different picture. Only 
additional observations, and longer time coverage, can tell us 
what is the case. With this caveat, this database is nevertheless a 
good starting point for the discussion of the TAC phenomenon 
in UX Ori. 

Among the detected events, the infalling gas (RACs) shows 
the largest velocities. For instance, event #1 shows a maximum 
average radial velocity of about 200 km/s which is roughly a 
half of the value of the stellar escape velocity, ~4 10 km/s, while 
the maximum average velocity of the outflowing gas (BACs) 
is about 100 km/s. This observational result could be biased 
since we begin detecting outflowing gas at maximum veloc- 
ities, which is not always the case for the infalling gas (e.g. 
events #4 and #7). We also note that blueshifted gas is detected 
when there are redshifted absorptions, but the opposite is not 
true; we do not know if there is a physical meaning behind this 
result. 

The dynamics of events #1 and #2 (RACs) and #3 and #5 
(BACs) denotes a deceleration of the gas. In event #4 the gas 
first accelerates and then decelerates; the data of #7 suggest an 
acceleration, while event #6 is only present in one spectrum. An 
estimate of the acceleration of the events over the time interval 
covered by the data is given in Table|. Positive values mean ac- 
celeration while negative ones mean deceleration. Table [| also 
gives the duration of each event, which is defined as the time it 
takes to reach zero velocity (for decelerating gas) or to go from 
zero to the maximum observed velocity for accelerating gas. 
The estimates suggest that the events typically last for a few 
days and that accelerations are a fraction of m s~ 2 . It is interest- 
ing to point out that the acceleration and deceleration phases of 
event #4 are of the same order of magnitude. We also note that 
our data do not allow us to conclude that this event is unique, 
since we may have missed the accelerating part of the other 
events. In this respect, it would have been of great interest to 
follow the evolution of event #7, which represents infalling gas 
with a significant acceleration, departing from a velocity very 
close to the stellar velocity (though it is only well detected in 
H/3 and Hy). 
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Fig. 5. Observed transient absorption components in each of the UES spectra (characterized by the Julian date of the observa- 
tions). Each point corresponds to one "TAC" and represents the mean velocity of the features with similar radial velocities de- 
tected in different lines. The radial velocity behaviour allows us to follow the temporal evolution of the gas causing the "TAC"". 
7 different dynamic events are identified; each event is enclosed in a box with an identification number (#1, #2,....). Error bars 
show the rms error of the average velocity and the number of lines used to estimate the average velocity is indicated. Fractional 
numbers come from the weight atributed to individual lines (see text). For the sake of clarity, the spectra taken during January 
31 st 1999 (JD 2451210.4) are expanded along the x-axis. [See the electronic edition for a colour version of this figure] 
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Table 3. Estimates of the acceleration and duration of each 
transient event. The time scale for the event's lifetime is es- 
timated as t = |V max /a|. The acceleration and deceleration 
phases of event #4 are shown separately, r for event #7 means 
the time needed to achieve V max departing from zero velocity. 



Event 


Type 


a(ms 2 ) 


t (days) 


1 


infall 


-0.8 


3.0 


2 


infall 


-0.1 


8.1 


3 


outflow 


-0.5 


2.1 


4+ 


infall 


+0.8 


2.6 


4- 


infall 


-0.8 


2.7 


5 


outflow 


-0.2 


3.2 


7 


infall 


+0.6 


0.3 



Table 4. Ratios of the average R parameter of several lines 
to Hy. The ratios (column 2) are estimated using all TACs in 
which the corresponding line absorption is present. A sigma- 
clipping algorithm has been applied in order to reject bad R 
values. The statistical error, the sigma-clipping value adopted 
and the fraction of the rejected values are given in columns 3 
to 5. The theoretical value for the Balmer lines is shown for 
comparison (gfii ne /gfH r )• 



Line 


Ratio 


Error 


clip, a 


% Rej. 


Theor. 


w 


1.2 


0.2 


2.5 


13 


2.67 


H<5 


0.8 


0.2 


2.5 





0.49 


H? 


0.6 


0.15 


2.5 





0.18 


CanK 


0.9 


0.2 


2.5 


5 




Nai D2 


0.14 


0.06 


2.0 


20 





All the transient line absorption features are very broad 
(Table |J), as is evident in the raw spectra. We have computed 
for each TAC of each event the average value of Av, weighted 
as in the case of v (the rms dispersion of Av among the dif- 
ferent lines in a TAC is typically ~25%). There is a tendency 
for Nai lines to be narrower than hydrogen and Can lines 
(Av NaI =s 0.5 AvH,can)- This will not be investigated further 
in this paper, but it could give interesting clues to the physical 
conditions of the gas. Fig. ^ shows Av versus v and Av/v versus 
v for each detected TAC, which are two complementary views. 
Each event with the corresponding TACs are enclosed in boxes. 
The diagrams show that events are well separated and that the 
velocity dispersion does not change drastically along their tem- 
poral evolution (perhaps with the exception of event #2). This 
behaviour is remarkable in event #4 which has both accelera- 
tion and deceleration phases. Thus, these results suggest that 
events could be characterized by a kind of "intrinsic" velocity 
dispersion throughout their lifetimes. 

4.2. R parameter 

Table || gives the parameter R for each absorption component. 
The ratio of the R values among different lines hardly varies 
for all 24 identified TACs. This is shown in Table [ij where 
we give the ratio of HyS, H6, Hf, Can K and Nai D2 to Hy. 
Hy has been chosen as the reference line because it is clearly 
present in all events. To estimate the errors we have applied a 
sigma-clipping algorithm to reject bad points; the correspond- 
ing values are given in Table || The remaining lines do not have 
a well defined ratio to Hy; we think this has no physical mean- 
ing, since He and Ca n H are blended and the relatively weak 
Na i Dl is more affected than Na i D2 by the strong telluric ab- 
sorptions in this wavelength interval. 

The fact that the ratio of R between different lines is rather 
constant allows us to characterize each TAC by an "intensity" 
(< R >) in the following manner. The ratios R\i n JRHy of each 
TAC are used to compute an "equivalent average Ru y < R >"■ 
< R > is estimated from Hfi, H6, Ca n K (weight 1 .0) and H£ 
(weight 0.5). Fig. || shows < R > as a function of the TAC 
velocities and events are enclosed in boxes. No correlation is 
seen. However, as in the case of Av, the estimated < R > values 



do not seem to change drastically from TAC to TAC within 
individual events. 



5. Discussion 

The intensity and dynamics of the transient components are re- 
vealed by the spectra in such detail that a careful analysis in 
the context of specific models (magnetospheric accretion mod- 
els, to mention one obvious example) could be used to validate 
the theoretical scenarios. This is in itself a huge effort, well be- 
yond the scope of this paper. There are, however, several points 
that we think are worth mentioning. Firstly, the data indicate 
that the events are the signatures of the dynamical evolution of 
gaseous clumps in the UX Ori CS disk. There do not seem to 
be large changes in the velocity dispersion and the relative line 
absorption strengths through the lifetime of each clump. This 
suggest that they are 'blobs' which basically preserve their ge- 
ometrical and physical identity. In this respect, event #4 is re- 
markable since during its estimated life of ~6 days it shows 
both an acceleration and a deceleration. In addition, the data 
reveal other aspects concerning the nature and dynamics of the 
events which are indicated in the following subsections. 

5.1. Origin of the circumstellar gas 

Based on spectra q ualitatively similar to those discussed here, 
Natta et aT| ( 2000| ) analysed the chemical composition of one 
infalling event in UX Ori. Non-LTE models were used to esti- 
mate the ionization and excitation of Balmer (up to H6), Ca II 
and Nai lines. Roughly solar abundances were found for the gas 
causing the redshifted features of that event. The present data 
suggest a similar nature for the gas causing the identified blobs. 
We detect evidence of infalling gas in the Balmer lines and in a 
number of metallic lines with similar velocities, which also dis- 
play roughly constant strength ratios to the Balmer lines; these 
ratios are < 1 (see Table Q). Similar result are found for the out- 
flowing gas (blueshifted absorptions). Thus, we can conclude 
that the CS gas in UX Ori is not very metal rich. In addition, 
the data suggest that the physical conditions of the outflowing 
and infalling gas are rather similar and they likely co-exist in 
space. 
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A second point is that in both RACs and BACs the CS gas 
has a significant underlying emission, at least in the Balmer 
lines we have examined. This can be seen easily if we compare 
the values of the observed average ratio for the Balmer lines 
to the ratio of the opacity in the lines (see Table |], columns 
2 and 6). If the CS gas was just absorbing the photospheric 
flux, then the ratio of any two lines originating from the same 
lower level should be equal to the ratio of their opacities. This 



is clearly not the case (see also similar results in Natta et al 



(2000)). There are a number of reasons why this can happen. 
The first is that the observed absorption is due to a very opti- 
cally thick cloud, whose projected size is smaller than the stel- 
lar surface. In this case, however, R should be the same for 
all the lines, roughly equal to the occulted fraction of the stel- 
lar surface, and this does not seem to be the case. In particu- 
lar, the Nai Dl absorption is always weaker than the Nai D2 
(RnsiD2 /Rnbidi = 1.4 + 0.4). A second possibility is that there 
is significant emission in the lines. 

R, for any given component, can be estimated from the ra- 
tio F b s /F syn . This can be done following eq. 5 of Rodgers et al 



( 2002 ), which provides this ratio under the assumption of a stel- 



lar and a circumstellar contribution to the observed flux; the cir- 
cumstellar contribution is caused by a spherical occulting cloud 
and a more extended envelope. In our zero order approach we 
neglect the extended envelope, and we assume a small opti- 
cal depth and black-body emission for the spherical occulting 
cloud. In this case, R can be written as: 



R 



B v (T ex ) //?cloud\ 2 

~ BAT,,) ' RZ ' 



(1) 



where t is the optical depth in the line, and T ex is the 
black-body excitation temperature of the gas in the cloud. R c \ ou d 
is the projected radius of the CS cloud. Some "filling-in" of 
the absorption features (R < r) may occur if r ex ~ and 
#cioud ~ R* (e.g. T ex a 7000/f and R c i ud/R* = 16, which is 
the corotation radius for UX Ori, reproduce the observed ra- 
tios). 

5.2. Gas Dynamics 

As already pointed out, the discussion in the previous subsec- 
tion together with the results of Natta et al. ( 2000 ) indicates that 



the TACs in UX Ori do not arise from the evaporation of solid 
bodies, but in a substantially different scenario. The simultane- 
ous presence of infalling and outflowing gas in the immediate 
vicinity of UX Ori is reminiscent of the predictions of magneto- 
spheric accretion models, as developed by a number of authors 



in the last several years (see Hartmann (1998) for a basic ac 



count of this theory). All these models assume stationary con- 
ditions, and it is difficult to derive quantitative predictions on 
the time dependent behaviour of the CS lines. However, some 
simple points can be made to constrain the different models. 

Let us start with the infall motions we detect. The sim- 
plest case we can consider for comparison is that of a blob 
of gas moving along the field lines of a stellar dipole mag- 
netic field, seen in absorption against the stellar photosphere. 
The infalling gas is practically contained in a region of the 



size of the corotation radius, which for UX Ori is R co ~ 1 .6 
R+. This is in agreement with the discussion in subsection 
5.1. The small value of the corotation radius implies that the 
maximum poloidal velocity the blob will reach is of order 
~ V esc (l - iW^co) ~ 220 kms" 1 . Since the UX Ori disk is 



seen almost edge-on ( Voshchinnikov et al. 1988), the observed 
velocity shift will be of the same order, which is in agreement 
with what we observe. A simple calculation shows that the blob 
will be seen to accelerate from very low velocity to this max- 
imum value in a time interval of about 1 day, until it hits the 
star. The acceleration seen in event #4 is reminiscent of such a 
behaviour, although the time scale is longer. On the contrary, 
the decelerating infall motions we see (events #1,2,4) require 
a modification of this simple model. A significant distortion of 
the velocity pattern can be caused by the rapid rotation of the 
star, which can create a sort of centrifugal barrier for the in- 
falling gas. Usin g the expression o f the poloidal velocity given 
by eq. (8.76) of ^artmann ( 1998 ) for the UX Ori parameters, 
we find a decrease of the maximum velocity (to about 180 
km s _1 ), but a similar acceleration pattern and time scale. Still, 
the interaction of the fast rotating magnetosphere with the ac- 
creting disk material should be considered in detail. Muzerolle 
et al. (2001 ) have included rotation in a study of T Tauri Stars 
(TTS), where, however, it is much smaller than the escape ve- 
locity at the stellar surface and the change of the poloidal veloc- 
ity is negligible. Similar calculations for cases where the stellar 
rotation is a significant fraction of the escape velocity would be 
very valuable. 

Outflowing gas is seen in events #3 and #5. In both cases 
we start detecting the outflow at negative velocities of 100- 
80 kms -1 , and we see deceleration on time scale of 2-3 days. 
Stellar-field-driven wind models (X-wind, see Shu et al. 2000, 



and references therein) predict that the wind leaves the disk sur- 
face close to the corotation radius on almost radial trajectories, 
accelerating rapidly to a speed of the order of Q+/? co ~ 350 
kms~'. This is significantly more than the maximum outflow 
velocity we observe, although the time coverage of our spectra 
can bias this result since we always find the maximum velocity 
of the outflow in the first spectrum taken in both the October 
98 and January 99 observing runs. Nevertheles, the decelera- 
tion we observe is not predicted by the models, unless colli- 
mation of the wind into the direction of the magnetic axis (i.e., 
away from the line of sight to the observer) occurs before the 
maximum speed is achieved and takes place on time scales of 
few day s. Wi nds driven by disk magnetic field (see Konigl & 
Pudritz 2000, and references therein) tend to produce slower 



outflows with large velocity di spersions from a broad region of 
the disk. |Goodson et al. ( 1997 ) discuss models where the inter- 
action between the rotating stellar field and a disk field results 
in a time-dependent launching of a two component outflow, of 
which one (the disk wind) has some similarity to what we ob- 
serve. However, also in these models it is not clear what could 
cause the observed deceleration. It is noteworthy that, once the 
acceleration phase is over, a wind can be decelerated by the ef- 
fect of the stellar gravity (flVlitskevich et al.| [l993|). In UX Ori 



-2 



the deceleration caused by the central star is ~ 0.2 - 0.5 ms 
at distances of ~ 14 - 20 R*, where the local escape velocity 
is ~ 130 kms -1 . An episode of outflow consistent with being 



10 



Mora et al.: A dynamical study of the circumstellar gas in UX Orionis 



gravitationally decelerated was observed in the TTS SU Aur by 
Petrov et al l (11996b. 



The predictions of simple magnetospheric accretion mod- 
els have been compared to the observed line variability in a 
number of TTS, with moderate success. The best case is that of 
SU Aur, a rather massive TTS with a rotation period of about 
3 days. This star shows evidence of simultaneous infall and 
outflow motions in various lines, similar to those observed in 



UX Ori (Johns & Basri 1995; Petrov et al. 1996; Oliveira et al 



200C). The infalling velocity is clearly modulated by the rota- 
tion period of the star, and a model with magnetically chanelled 
accretion in a dipole field inclined by a few degrees to the ro- 
tation axis of the star may account for most of the observed 
infall properties. As in UX Ori the outflow is very variable, 
with maximum velocity similar to that of the infalling gas and 
a tendency to decelerate o n a time scale of days , with no obvi- 
ous rotational modulation ( Johns & Basri 1995 ). Note that the 
timing strategy of our UX Ori observations is not suited to the 
detection of modulations on the time scale of the stellar rotation 
period, which is about 17 hr. 

UX Ori may provide a useful test-case. On one hand, it is 
likely that the structure of the stellar magnetic field is more 
complex than in TTS, since UX Ori (a A-type star of 2-3 M ) 
lacks the surface convective layer that may create a simple 
magnetic field in TTS. On the other hand, one may take ad- 
vantage of the fact that the average accretion rate of UX Ori is 
rather low ( Tambovtseva et al. 2001 ). The CS gas is mostly seen 
in absorption against the continuum, when a sporadic increase 
of the accretion rate in the disk creates gas "blobs" which move 
along the magnetic field lines, tracing their pattern more clearly 
than in objects with higher accretion rates, where broad emis- 
sion dominates the CS lines. We would like to mention that 



Pontefract et al. (2000) also suggest magnetospheric accretion 



as a promising model to explain the Ha spectropolarimetry data 
of the Herbig Ae star AB Aur. 

One final point to keep in mind in this context is that varia- 
tions along the flow (both for infall and outflow) of the source 
function in the various lines may have an important effect on 
the comparison between observations and models, so that de- 
tailed calculations are required for validation of the models. 

6. Conclusions 

The data presented in this paper allow us to analyse the spec- 
troscopic behaviour of the CS gas disk around UX Ori on time 
scales of months, days and hours. Significant activity in the CS 
disk is always present, which manifests itself in the continuous 
appearance and disappearance of absorption components de- 
tected in hydrogen and in many metallic lines. This activity is 
not related to substantial variations of the stellar photosphere. 
Blobs of gas experiencing infalling and outflowing motions are 
the likely origin of the transient features. Blobs undergo accel- 
erations/decelerations of the order of tenths of m s~ 2 and last 
for a few days. Detectable changes in the gas dynamics occur 
on a time scale of hours, but the intrinsic velocity dispersion 
of the blobs appears to remain rather constant. No noticeable 
differences are seen in the properties of the infalling and out- 
flowing gas, although infalls might have larger velocity disper- 



sion. The relative absorption strength of the transient absorp- 
tions suggests gas abundances similar to the solar metallicity, 
ruling out the evaporation of solid bodies as the physical origin 
of the spectroscopic features. We suggest that the data should 
be analysed in the context of detailed magnetospheric accretion 
models, similar to those used for T Tauri stars. 
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